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QQ ' Abstract 

"^ ! The dynamical generation of a four- dimensional classical universe from nothing 

^ I but fundamental quantum excitations at the Planck scale is a long-standing chal- 

(yQ I lenge to theoretical physicists. A candidate theory of quantum gravity which 

O ' achieves this goal without invoking exotic ingredients or excessive fine-tuning is 

>• ■ based on the nonperturbative and background-independent technique of Causal 

K> , Dynamical Triangulations. We demonstrate in detail how in this approach a 

H I macroscopic de Sitter universe, accompanied by small quantum fluctuations, 

■ ■ ■ ' emerges from the full gravitational path integral, and how the effective action 

determining its dynamics can be reconstructed uniquely from Monte Carlo data. 

We also provide evidence that it may be possible to penetrate to the sub-Planckian 

regime, where the Planck length is large compared to the lattice spacing of the 

underlying regularization of geometry. 



1 Introduction 

A major unsolved problem in theoretical physics is to reconcile the classical the- 
ory of general relativity with quantum mechanics. Of the numerous attempts, 
some have postulated new and so far unobserved ingredients, while others have 
proposed radically new principles governing physics at the as yet untested Planck- 
ian energy scale. Here we report on a much more mundane approach using only 
standard quantum field theory. In a sum-over-histories approach we will attempt 
to define a nonperturbative quantum field theory which has as its infrared limit 
ordinary classical general relativity and at the same time has a nontrivial ultra- 
violet limit. From this point of view it is in the spirit of the renormalization 
group approach, first advocated long ago by Weinberg [1], and more recently 
substantiated by several groups of researchers [2] . However, it has some advan- 
tages compared to the renormalization group approach in that it allows us to 
study (numerically) certain geometric observables which are difficult to handle 
analytically. 

We define the path integral of quantum gravity nonperturbatively using the 
lattice approach known as causal dynamical triangulations (CDT) as a regular- 
ization. In Sec. 2 we give a short description of the formalism, providing the 
definitions which are needed later to describe the measurements. CDT estab- 
lishes a nonperturbative way of performing the sum over four-geometries (for 
more extensive definitions, see [3, 4]). It sums over the class of piecewise linear 
four-geometries which can be assembled from four-dimensional simplicial building 
blocks of link length a, such that only causal spacetime histories are included. 
The continuum limit of such a lattice theory should ideally be obtained as for 
QCD defined on an ordinary fixed lattice, where for an observable (9(x„), Xn 
denoting a lattice point, one can measure the correlation length ^{go) from 

- \og{{0{xr,)0{ym))) ~ |n - m|/e((?o) + o{\n-m\). (1) 

A continuum limit of the lattice theory may then exist if it is possible to fine- 
tune the bare coupling constant qq of the theory to a critical value g^ such that 
the correlation length goes to infinity, ^{go) -^ oo. Knowing how ^{go) diverges 
for go -^ g^ determines how the lattice spacing a should be taken to zero as a 
function of the coupling constants, namely 

^(^o) = ] ^, a{go) = \go - gX- (2) 

\9o — 9{)\ 

The challenge when searching for a field theory of quantum gravity is to find 
a theory which behaves in this way. The challenge is three-fold: (i) to find a 
suitable nonperturbative formulation of such a theory which satisfies a minimum 
of reasonable requirements, (ii) to find observables which can be used to test 



relations like (1), and (iii) to show that one can adjust the coupling constants 
of the theory such that (2) is satisfied. Although we will focus on (i) in what 
follows, let us immediately mention that (ii) is notoriously difficult in a theory of 
quantum gravity, where one is faced with a number of questions originating in the 
dynamical nature of geometry. What is the meaning of distance when integrating 
over all geometries? How do we attach a meaning to local spacetime points like 
Xn and |/„? How can we define at all local, diffeomorphism-invariant quantities in 
the continuum which can then be translated to the regularized (lattice) theory? 
- What we want to point out here is that although (i)-(iii) are standard require- 
ments when relating critical phenomena and (Euclidean) quantum field theory, 
gravity is special and may require a reformulation of (part of) the standard sce- 
nario sketched above. We will return to this issue when we discuss our results in 
Sec. 8. 

Our proposed nonperturbative formulation of four-dimensional quantum grav- 
ity has a number of nice features. Firstly, it sums over a class of piecewise linear 
geometries, which - as usual - are described without the use of coordinate sys- 
tems. In this way we perform the sum over geometries directly, avoiding the 
cumbersome procedure of first introducing a coordinate system and then getting 
rid of the ensuing gauge redundancy, as one has to do in a continuum calculation. 
Our underlying assumptions are that 1) the class of piecewise linear geometries is 
in a suitable sense dense in the set of all geometries relevant for the path integral 
(probably a fairly mild assumption), and 2) that we are using a correct measure 
on the set of geometries. This is a more questionable assumption since we do 
not even know whether such a measure exists. Here one has to take a pragmatic 
attitude in order to make progress. We will simply examine the outcome of our 
construction and try to judge whether it is promising. 

Secondly, our scheme is background-independent. No distinguished geometry, 
accompanied by quantum fiuctuations, is put in by hand. If the CDT-regularized 
theory is to be taken seriously as a potential theory of quantum gravity, there 
has to be a region in the space spanned by the bare coupling constants where 
the geometry of spacetime bears some resemblance with the kind of universe 
we observe around us. That is, the theory should create dynamically an effective 
background geometry around which there are (small) quantum fluctuations. This 
is a very nontrivial property of the theory and one we are going to investigate 
in detail in the present piece of work. New computer simulations presented here 
confirm in a much more direct way the indirect evidence for such a scenario which 
we provided earlier in [6, 7]. They establish the de Sitter nature of the background 
spacetime, quantify the fiuctuations around it, and set a physical scale for the 
universes we are dealing with. The main results of our investigation, without the 
numerical details, were announced in [8] (see also [9]). 

The rest of the article is organized as follows. In Sec. 2 we describe briefiy 



the regularization method of quantum gravity named CDT and the set-up of the 
computer simulations. In Sec. 3 we present the evidence for an effective back- 
ground geometry corresponding to the four-dimensional sphere S^, i.e. Euclidean 
de Sitter spacetime. Sec. 4 deals with the reconstruction of an effective action 
for the scale factor of the universe from the computer data, and in Sec. 5 we 
analyze the quantum fluctuations around the "classical" S'^-solution. Sec. 6 con- 
tains an analysis of the geometry of the spatial slices of our computer-generated 
universe. In Sec. 7 we determine the physical sizes of our universes expressed in 
Planck lengths and try to follow the flow of the gravitational coupling constant 
of the effective action under a change of the bare coupling constants of the bare, 
"classical" action used in the path integral. Finally we discuss the results, their 
interpretation and future perspectives of the CDT-quantum gravity theory in Sec. 
8. 

2 Causal Dynamical Triangulations (CDT) 

The approach of causal dynamical triangulations stands in the tradition of [11], 
which advocated that in a gravitational path integral with the correct, Lorentzian 
signature of spacetime one should sum over causal geometries only. More specif- 
ically, we adopted this idea when it became clear that attempts to formulate a 
Euclidean nonperturbative quantum gravity theory run into trouble in spacetime 
dimension d larger than two. At the same time, such a causal reformulation re- 
sults in a path integral which relates more closely to canonical formulations of 
quantum gravity. 

This implies that we start from Lorentzian simplicial spacetimes with d = 4 
and insist that only causally well-behaved geometries appear in the (regularized) 
Lorentzian path integral. A crucial property of our explicit construction is that 
each of the configurations allows for a rotation to Euclidean signature. We rotate 
to a Euclidean regime in order to perform the sum over geometries (and rotate 
back again afterwards if needed). We stress here that although the sum is per- 
formed over geometries with Euclidean signature, it is different from what one 
would obtain in a theory of quantum gravity based ab initio on Euclidean space- 
times. The reason is that not all Euclidean geometries with a given topology are 
included in the "causal" sum since in general they have no correspondence to a 
causal Lorentzian geometry. 

How do we construct the class of piecewise linear geometries used in the 
Lorentzian path integral (see [3] for a detailed description)? The most important 
assumption is the existence of a global proper-time foliation. We assume that 
the spacetime topology is that of / x T.^^\ where S*^^^ denotes an arbitrary three- 
dimensional manifold. In what follows, we will for simplicity study the case of 



the simplest spatial topology T.^^^ = S^, that of a three-sphere. The compactness 
of S^ obviates the discussion of spatial boundary conditions for the universe. 
The spatial geometry at each discrete proper-time step t„ is represented by a 
triangulation of S^, made up of equilateral spatial tetrahedra with squared side- 
length £^ = a^ > 0. In general, the number A^3(t„) of tetrahedra and how 
they are glued together to form a piecewise flat three-dimensional manifold will 
vary with each time-step t„. In order to obtain a four- dimensional triangulation, 
the individual three-dimensional slices must still be connected in a causal way, 
preserving the S'^-topology at all intermediate times t between t„ and tn+i-^ This 
is done by connecting each tetrahedron belonging to the triangulation at time t„ 
to a vertex belonging to the triangulation at time tn+i by means of a four-simplex 
which has four time-like links of length-squared if = —cxd-l-, a > 0, interpolating 
between the adjacent slices (a so-called (4,l)-simplex). In addition, a triangle in 
the triangulation at time t„ can be connected to a link in the triangulation at 
tn+i via a four-simplex with six time-like links (a so-called (3,2)-simplex), again 
with dl = —aif. Conversely, one can connect a link at t„ to a triangle at tn+i 
to create a (2,3)-simplex and a vertex at t„ to a tetrahedron at t„+i to create a 
(l,4)-simplex. One can interpolate between subsequent triangulations of S^ at 
tn and tn+i in many distinct ways compatible with the topology J x S*^ of the 
four-manifold. All these possibilities are summed over in the CDT path integral. 
The explicit rotation to Euclidean signature is done by performing the rotation 
a -^ —a in the complex lower half-plane, \a\ > 7/12, such that we have if = \a\if 
(see [3] for a discussion). 

The Einstein-Hilbert action S^^ has a natural geometric implementation on 
piecewise linear geometries in the form of the Regge action. This is given by the 
sum of the so-called deficit angles around the two-dimensional "hinges" (subsim- 
plices in the form of triangles), each multiplied with the volume of the correspond- 
ing hinge. In view of the fact that we are dealing with piecewise linear, and not 
smooth metrics, there is no unique "approximation" to the usual Einstein-Hilbert 
action, and one could in principle work with a different form of the gravitational 
action. We will stick with the Regge action, which takes on a very simple form 
in our case, where the piecewise linear manifold is constructed from just two dif- 
ferent types of building blocks. After rotation to Euclidean signature one obtains 



"'^This implies the absence of branching of the spatial universe into several disconnected 
pieces, so-called baby universes, which (in Lorentzian signature) would inevitably be associated 
with causality violations in the form of degeneracies in the light cone structure, as has been 
discussed elsewhere (see, for example, [10]). 



for the action (see [4] for details) 



^^Regge ^ _(^^ + 6A)iVo + /s:4(iVf'^+iVf'^) + A(2iVf')+iVf')), (3) 

where A^o denotes the total number of vertices in the four-dimensional triangula- 
tion and AQ ' and AQ ' denote the total number of the four-simplices described 
above, i.e. the total number of (4,l)-simplices plus (l,4)-simplices and the to- 
tal number of (3,2)-simplices plus (2,3)-simplices, respectively, so that the total 
number N4 of four-simplices is A^4 = A'4 + A^] . The dimensionless coupling 
constants kq and ^4 are related to the bare gravitational and bare cosmologi- 
cal coupling constants, with appropriate powers of the lattice spacing a already 
absorbed into Kq and k,^. The asymmetry parameter A is related to the parame- 
ter a introduced above, which describes the relative scale between the (squared) 
lengths of space- and time-like links. It is both convenient and natural to keep 
track of this parameter in our set-up, which from the outset is not isotropic in 
time and space directions, see again [4] for a detailed discussion. Since we will 
in the following work with the path integral after Wick rotation, let us redefine 
a := —a [4], which is positive in the Euclidean domain.^ For future reference, 
the Euclidean four-volume of our universe for a given choice of a is given by 

4 ^V8a-3 (4,1) Vl2a - 7 (3,2) 



V, = C, a' (^^^^ Nr' + ^ ^ Nr') , (4) 

where C4 = v5/96 is the four-volume of an equilateral four-simplex with edge 
length a = 1 (see [3] for details). It is convenient to rewrite expression (4) as 

V, = C,{0 «' ATf '^ = ^4(0 a' N,/{1 + 0, (5) 

where ^ is the ratio 

C = Nr^/Nt\ (6) 

and (74(^)0^ is a measure of the "effective four- volume" of an "average" four- 
simplex. In computing (3), we have assumed that the spacetime manifold is com- 
pact without boundaries, otherwise appropriate boundary terms must be added 
to the action. 

The path integral or partition function for the CDT version of quantum gravity 
is now 

Z(G,A) = y"l?Me-^i"[^l ^ Z(/.o,«:4,A) = 5]^e-^-(^), (7) 



^The most symmetric choice is a = 1, corresponding to vanishing asymmetry, A = 0. 



where the summation is over all causal triangulations T of the kind described 
above, and we have dropped the superscript "Regge" on the discretized action. 
The factor 1/Cr is a symmetry factor, given by the order of the automorphism 
group of the triangulation T. The actual set-up for the simulations is as follows. 
We choose a fixed number A^ of spatial slices at proper times ti, t2 = ti + a^, up 
to In = ti + {N — l)at., where At = at is the discrete lattice spacing in temporal 
direction and T = Nat the total extension of the universe in proper time. For 
convenience we identify tN+i with ti, in this way imposing the topology S^ x S^ 
rather than I x S^. This choice does not affect physical results, as will become 
clear in due course. 

Our next task is to evaluate the nonperturbative sum in (7), if possible, an- 
alytically. Although this can be done in spacetime dimension d = 2 ([5], and 
see [12] for recent developments) and at least partially in d = 3 [13, 14], an an- 
alytic solution in four dimensions is currently out of reach. However, we are in 
the fortunate situation that Z{k,q, k,4, A) can be studied quantitatively with the 
help of Monte Carlo simulations. The type of algorithm needed to update the 
piecewise linear geometries has been around for a while, starting from the use of 
dynamical triangulations in bosonic string theory (two-dimensional Euclidean tri- 
angulations) [15, 16, 17] and was later extended to their application in Euclidean 
four-dimensional quantum gravity [18, 19]. In [3] the algorithm was modified to 
accomodate the geometries of the CDT set-up. Note that the algorithm is such 
that it takes the symmetry factor Cr into account automatically. 

We have performed extensive Monte Carlo simulations of the partition func- 
tion Z for a number of values of the bare coupling constants. As reported in [4] , 
there are regions of the coupling constant space which do not appear relevant for 
continuum physics in that they seem to suffer from problems similar to the ones 
found earlier in Euclidean quantum gravity constructed in terms of dynamical tri- 
angulations, which essentially led to its abandonement in d > 2. Namely, when 
the (inverse, bare) gravitational coupling k,q is sufficiently large, the Monte Carlo 
simulations exhibit a sequence in time direction of small, disconnected universes, 
none of them showing any sign of the scaling one would expect from a macro- 
scopic universe. We believe that this phase of the system is a Lorentzian version 
of the branched polymer phase of Euclidean quantum gravity. By contrast, when 
A is sufficiently small the simulations reveal a universe with a vanishing tempo- 
ral extension of only a few lattice spacings, ending both in past and future in a 
vertex of very high order, connected to a large fraction of all vertices. This phase 
is most likely related to the so-called crumpled phase of Euclidean quantum grav- 
ity. The crucial and new feature of the quantum superposition in terms of causal 
dynamical triangulations is the appearance of a region in coupling constant space 
which is different and interesting and where continuum physics may emerge. It 
is in this region that we have performed the simulations reported in this article. 



and where previous work has aheady uncovered a number of intriguing physical 
results [6, 7, 4, 20]. 

In the Euclideanized setting the value of the cosmological constant determines 
the spacetime volume V4 since the two appear in the action as conjugate vari- 
ables. We therefore have (V4) ~ G/A in a continuum notation, where G is the 
gravitational coupling constant and A the cosmological constant. In the com- 
puter simulations it is more convenient to keep the four-volume fixed or partially 
fixed. We will implement this by fixing the total number of four-simplices of type 
AQ ' or, equivalently, the total number N^^ of tetrahedra making up the spatial 
S^ triangulations at times tj, i = 1, . . . , A^, 



N 1 

N, = Y.N,{U) = -N^t"^. (8) 



We know from the simulations that in the phase of interest (AQ ) oc (AQ ) as 
the total volume is varied [4]. This effectively implies that we only have two bare 
coupling constants kq, A in (7), while we compensate by hand for the coupling 
constant K/^ by studying the partition function Z{kq, A; AQ ' ) for various A^, ' 



4 



To keep track of the ratio ^(kq, A) between the expectation value (AQ ' ) and 
A4 , which depends weakly on the coupling constants, we write (c.f. eq. (6)) 

(N,) = ATf ^) + (A^f ')) = ATf ^)(1 + e(«:o, A)). (9) 

For all practical purposes we can regard A4 in a Monte Carlo simulation as fixed. 
The relation between the partition function we use and the partition function 
with variable four-volume is given by the Laplace transformation 

POD 

Z{Ko,Ki,A)= dNie-^"^" Z{Ko,Ni,A), (10) 

Jo 

where strictly speaking the integration over A4 should be replaced by a summa- 
tion over the discrete values A4 can take. 



3 The macroscopic de Sitter universe 

The Monte Carlo simulations referred to above will generate a sequence of space- 
time histories. An individual spacetime history is not an observable, in the same 
way as a path x{t) of a particle in the quantum-mechanical path integral is not. 
However, it is perfectly legitimate to talk about the expectation value {x{t)) as 
well as the fluctuations around {x{t)). Both of these quantities are in principle 
calculable in quantum mechanics. 



Obviously, there are many more dynamical variables in quantum gravity than 
there are in the particle case. We can still imitate the quantum-mechanical situ- 
ation by picking out a particular one, for example, the spatial three- volume Vs{t) 
at proper time t. We can measure both its expectation value (V3(t)) as well as 
fluctuations around it. The former gives us information about the large-scale 
"shape" of the universe we have created in the computer. In this section, we will 
describe the measurements of (V3(t)), keeping a more detailed discussion of the 
fluctuations to Sec. 5 below. 

A "measurement" of V3(t) consists of a table N^ii), where i = 1, . . . ,N denotes 
the number of time-slices. Recall from Sec. 2 that the sum over slices Yli=i ^3(0 
is kept constant. The time axis has a total length of N time steps, where A^ = 80 
in the actual simulations, and we have cyclically identifled time-slice N + 1 with 
time-slice 1. 

What we observe in the simulations is that for the range of discrete volumes 
A4 under study the universe does not extend (i.e. has appreciable three-volume) 
over the entire time axis, but rather is localized in a region much shorter than 
80 time slices. Outside this region the spatial extension N^li) will be minimal, 
consisting of the minimal number (flve) of tetrahedra needed to form a three- 
sphere S^, plus occasionally a few more tetrahedra.^ This thin "stalk" therefore 
carries little four-volume and in a given simulation we can for most practical 
purposes consider the total four- volume of the remainder, the extended universe, 
as flxed. 

In order to perform a meaningful average over geometries which explicitly 
refers to the extended part of the universe, we have to remove the translational 
zero mode which is present. During the Monte Carlo simulations the extended 
universe will fluctuate in shape and its centre of mass (or, more to the point, its 
centre of volume) will perform a slow random walk along the time axis. Since we 
are dealing with a circle (the compactifled time axis), the centre of volume is not 
uniquely deflned (it is clearly arbitrary for a constant volume distribution), and 
we must flrst deflne what we mean by such a concept. Here we take advantage of 
the empirical fact that our dynamically generated universes decompose into an 
extended piece and a stalk, with the latter containing less than one per cent of 
the total volume. We are clearly interested in a deflnition such that the centre of 
volume of a given conflguration lies in the centre of the extended region. One also 
expects that any sensible deflnition will be unique up to contributions related to 
the stalk and to the discreteness of the time steps. In total this amounts to an 
ambiguity of the centre of volume of one lattice step in the time direction. 

In analyzing the computer data we have chosen one speciflc deflnition which 



■^This kinematical constraint ensures that the triangulation remains a simplicial manifold in 
which, for example, two d-simpHces are not allowed to have more than one (d — l)-simplex in 



common. 



is in accordance with the discussion above^. Maybe surprisingly, it turns out 
that the inherent ambiguity in the choice of a definition of the centre of volume 
- even if it is only of the order of one lattice spacing - will play a role later 
on in our analysis of the quantum fiuctuations. For each universe used in the 
measurements (a "path" in the gravitational path integral) we will denote the 
centre-of-volume time coordinate calculated by our algorithm by icv From now 
on, when comparing different universes, i.e. when performing ensemble averages, 
we will redefine the temporal coordinates according to 

N^^^it) = N,{1 + mod(^ + 2e. - 1, N)), (11) 

such that the centre of volume is located at 0. 

Having defined in this manner the centre of volume along the time-direction of 
our spacetime configurations we can now perform superpositions of such config- 
urations and define the average {Ns{i)) as a function of the discrete time i. The 
results of measuring the average discrete spatial size of the universe at various 
discrete times i are illustrated in Fig. 1 and can be succinctly summarized by the 
formula 

Nii,) .- iNm - 5i^ ^_L_ ,„,. ^__ I . , . 0,59. ,12) 

where N^li) denotes the number of three-simp lices in the spatial slice at dis- 
cretized time i and A^4 the total number of four-simplices in the entire universe. 
Since we are keeping A^] ' fixed in the simulations and since ^ changes with the 
choice of bare coupling constants, it is sometimes convenient to rewrite (12) as 

mfii) = -Ni"''^ ^- cos^ I 1 , (13) 

where sq is defined by so(AQ ' )^/'' = sqN^ . Of course, formula (12) is only 
valid in the extended part of the universe where the spatial three-volumes are 
larger than the minimal cut-off size. 



^Explicitly, we consider the quantity 
CV{i') = 



N/2-1 

Y^ (i + 0.5) Nail + iaod{i' + i-l,N)) 

i=-N/2 



and find the value of i' e {1, . . . , iV} for which CV{i') is smallest. We denote this i' by icv 
If there is more than one minimum, we choose the value which has the largest three-volume 
N^ii'). Let us stress that this is just one of many definitions of 1^. All other sensible definitions 
will for the type of configurations considered here agree to within one lattice spacing. 
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Figure 1: Background geometry (N^^i)): MC measurements for fixed A'^^ = 
160.000 (A^4 = 362.000) and best fit (12) yield indistinguishable curves at given 
plot resolution. The bars indicate the average size of quantum fluctuations. 

The data shown in Fig. 1 have been collected at the particular values (kq, A) = 
(2.2, 0.6) of the bare coupling constants and for A^4 = 362.000 (corresponding to 
AQ = 160.000). For these values of (kq, A) we have verified relation (12) for 
N4 ranging from 45.500 to 362.000 building blocks (45.500, 91.000, 181.000 and 
362.000). After rescaling the time and volume variables by suitable powers of 
A^4 according to relation (12), and plotting them in the same way as in Fig. 
1, one finds almost total agreement between the curves for different spacetime 
volumes.^ Eq. (12) shows that spatial volumes scale according to A^4 and time 
intervals according to A^4 , as one would expect for a genuinely /ottr-dimensional 
spacetime. This strongly suggests a translation of (12) to a continuum notation. 
The most natural identification is given by 



WtV,^\t) 



\/4— COS^ 



where we have made the identifications 
ti i 



B 



soK 



1/4 ■ 



Ativ^ VsiU) = 2C^Ns{t)a\ 



such that we have 



dty^tV3{t)=Vi. 



(14) 



(15) 



(16) 



^By contrast, the quantum fluctuations indicated in Fig. 1 as vertical bars are volume- 
dependent and will be the larger the smaller the total four-volume, see Sec. 5 below for details. 



11 



In (15), ^Jgtt is the constant proportionality factor between the time t and genuine 
continuum proper time t^ t = ^Jgu t. (The combination Atjy^Va contains 
(74, related to the four-volume of a four-simplex rather than the three-volume 
corresponding to a tetrahedron, because its time integral must equal V4). Writing 
V4 = 8n'^R^/3, and y^ = R/B, eq. (14) is seen to describe a Euclidean de Sitter 
universe (a four-sphere, the maximally symmetric space for positive cosmological 
constant) as our searched-for, dynamically generated background geometry! In 
the parametrization of (14) this is the classical solution to the action 



-.tti'' ^' 



1 /■ .. ^ / 3'% (t) . , „,/3, 



'-2^ J *^ (^#f ^ *^'" *" - '''=<'V • "'' 

where k2 = 9(27r^)^/^ and A is a Lagrange multiplier, fixed by requiring that the to- 
tal four- volume be V4, / dt^/gu V^{t) = V4. Up to an overall sign, this is precisely 
the Einstein-Hilbert action for the scale factor a{t) of a homogeneous, isotropic 
universe (rewritten in terms of the spatial three-volume V^lt) = 27r^a(t)^), al- 
though we of course never put any such simplifying symmetry assumptions into 
the CDT model. 

For a fixed, finite four-volume V4 and when applying scaling arguments it 
can be convenient to rewrite (17) in terms of dimensionless units by introducing 
s = t/V^ and V3(t) = V^ Vs{s), in which case (17) becomes 

now assuming that f ds^/gZ v^^s) = 1, and with Qss = gtt- A discretized, dimen- 
sionless version of (17) is 

Sdiscr = ki2_^\ j;^. + hN^' {i) j , (19) 

where k2 oc A;2- This can be seen by applying the scaling (12), namely, Ns{i) = 
NJ n^{si) and Sj = i/Nj . With this scaling, the action (19) becomes 

S,.s.r = A:,v^$: As L^^ ^MfmW^s^y ^ A:,nf (s.)") , (20) 

where As = 1/N^/^^ and therefore has the same form as (18). This enables us 
to finally conclude that the identifications (15) when used in the action (19) lead 
naively to the continuum expression (17) under the identification 

^^a VG4^^ (21) 

ki 3^6 
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Figure 2: The measured average shape {N^{i)) of the quantum universe at A = 
0.6, for kq = 2.2 (broader distribution) and kq = 3.6 (narrower distribution) 



taken at N' 



(4,1) 



160.000. 



Next, let us comment on the universahty of these results. First, we have 
checked that they are not dependent on the particular definition of time-slicing 
we have been using, in the following sense. By construction of the piecewise 
linear CDT-geometries we have at each integer time step ti = iat a. spatial surface 
consisting of A^3(z) tetrahedra. Alternatively, one can choose as reference slices for 
the measurements of the spatial volume non-integer values of time, for example, 
all time slices at discrete times 2 — 1/2, i = 1,2, ... . In this case the "triangulation" 
of the spatial three-spheres consists of tetrahedra - from cutting a (4,1)- or a (1,4)- 
simplex half-way - and "boxes", obtained by cutting a (2,3)- or (3,2)-simplex 
(the geometry of this is worked out in [21]). We again find a relation like (12) 
if we use the total number of spatial building blocks in the intermediate slices 
(tetrahedra-(-boxes) instead of just the tetrahedra. 

Second, we have repeated the measurements for other values of the bare cou- 
pling constants. As long as we stay in the phase where an extended universe is 
observed (called "phase C" in ref. [4]), a relation like (12) remains valid. In addi- 
tion, the value of Sq, defined in eq. (12), is almost unchanged until we get close to 
the phase transition lines beyond which the extended universe disappears. Fig. 
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2 shows the average shape (A^3(t)) for A = 0.6 and for kq equal to 2.2 and 3.6. 
Only for the values of Kq around 3.6 and larger will the measured (N^lt)) differ 
significantly from the value at 2.2. For values larger than 3.8 (at A = 0.6), the 
universe will disintegrate into a number of small and disconnected components 
distributed randomly along the time axis, and one can no longer fit the distri- 
bution {N^it)) to the formula (12). Fig. 3 shows the average shape (N^^t)) for 
kq = 2.2 and A equal to 0.2 and 0.6. Here the value A = 0.2 is close to the 
phase transition where the extended universe will flatten out to a universe with 
a time extension of a few lattice spacings only. Later we will show that while Sq 
is almost unchanged, the constant ki in (19), which governs the quantum fluctu- 
ations around the mean value (N^lt)), is more sensitive to a change of the bare 
coupling constants, in particular in the case where we change kq (while leaving 
A fixed). 
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Figure 3: The measured average shape {N^^i)) of the quantum universe at kq = 
2.2, for A = 0.6 (broad distribution) and A= 0.2 (narrow distribution), both 



taken at N, 



(4,1) 



160.000. 
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4 Constructive evidence for the effective action 

While the functional form (12) for the three- volume fits the data perfectly and the 
corresponding continuum effective action (17) reproduces the continuum version 
(14) of (12), it is still of interest to check to what extent one can reconstruct 
the discretized version (19) of the continuum action (17) from the data explicitly. 
Stated differently, we would like to understand whether there are other effective 
actions which reproduce the data equally well. As we will demonstrate by explicit 
construction in this section, there is good evidence for the uniqueness of the action 
(19). 

The data we have are two-fold: the measurement of Ns{i), that is, the three- 
volume at the discrete time step i, and the measurement of the three-volume 
correlator Ns{i)Ns{j). Having created K statistically independent configurations 
AQ (i) by Monte Carlo simulation allows us to construct the average 

N,{^)■.= {N,{^)) = ^J2^^'^('^^ (22) 

k 

where the superscript in (■)*''^^ denotes the result of the k'th configuration sampled, 
as well as the covariance matrix 

^(^'•^■) = lT.(^t\^ - m)){Nt\j)-mj))- (23) 



k 



Since we have fixed the sum Yli=i ^3(0 (recall that A^ denotes the fixed number 
of time steps in a given simulation), the covariance matrix has a zero mode, 
namely, the constant vector e\ , 

J2C{t,j)ef = 0, ef^ = l/VN Vi (24) 

i 

A spectral decomposition of the symmetric covariance matrix gives 

C'=5^A,|e("))(e(«)|, (25) 



a=l 



where we assume the A^ — 1 other eigenvalues of the covariance matrix Cij are 
different from zero. We now define the "propagator" P as the inverse of C on the 
subspace orthogonal to the zero mode e^^^ , that is, 

N-l 

P=^^|e('^))(e('^)| = (C' + i)-i-i, i= |e(°))(e(°)|. (26) 



a=l 
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We now assume we have a discretized action which can be expanded around 
the expectation value Nsii) according to 

Sdzscr [N + n] = Sd^scr W] + 'Y, ^^^'^i + 0{n^). (27) 

If the quadratic approximation describes the quantum fluctuations around the 
expectation value N well, the inverse of P will be a good approximation to the 
covariance matrix. Conversely, still assuming the quadratic approximation gives 
a good description of the fluctuations, the P constructed from the covariance 
matrix will to a good approximation allow us to reconstruct the action via (27). 
Simply by looking at the inverse P of the measured covariance matrix, defined 
as described above, we observe that it is to a very good approximation small and 
constant except on the diagonal and the entries neighbouring the diagonal. We 
can then decompose it into a "kinetic" and a "potential" term. The kinetic part 
pkm jg defined as the matrix with non-zero elements on the diagonal and in the 
neighbouring entries, such that the sum of the elements in a row or a column is 
always zero, 

N 
pkin^J2p,X(i)^ (28) 

i=l 

where the matrix X^''' is given by 

-^jk ~ ^ij^ik + ^{i+l)j^{i+l)k — %+l)j'^Jfc — SijS{i+l)k- (29) 

Note that the range of P^*" lies by definition in the subspace orthogonal to the 
zero mode. Similarly, we define the potential term as the projection of a diagonal 
matrix D on the subspace orthogonal to the zero mode 

N 
ppot ^ (J _ J^^f)^f _ ^) _ ^ ^.y W. (30) 

j=l 

The diagonal matrix D and the matrices Y^"^' are defined by 

Djk = UjSjk, YJl^ = 5ij5ik - "^ ^ "^ + ^, (31) 

and / denotes the N x N unit matrix. 

The matrix P is obtained from the numerical data by inverting the covariance 
matrix C after subtracting the zero mode, as described above. We can now try 
to find the best values of the pj's and -Uj's by a least-x^ fit^ to 

tr (P - (P'^^" + PP°*)V . (32) 
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Figure 4: The directly measured expectation values Nsii) (thick gray curves), 
compared to the averages N3{i) reconstructed from the measured covariance ma- 
trix C (thin black curves), for kq = 2.2 and A = 0.6, at various fixed volumes 
N^ . The two-fold symmetry of the interpolated curves around the central 
symmetry axis results from an explicit symmetrization of the collected data. 

Let us look at the discretized minisuperspace action (19) which obviously has 
served as an inspiration for the definitions of P*^*" and pp°*. Expanding N^li) to 
second order around N3(i) one obtains the identifications 



mi) 



2ki 
Pi ' 



U"{N3i^)) = -u^ 



(33) 



rV3 



where f/(A^3(i)) = kik2N^ (i) denotes the potential term in (19). We use the 
fitted coefficients pi to reconstruct Nsi^i) and then compare these reconstructed 



^A x^"fit of the form (32) gives the same weight to each three- volume A^3(i). One might 
argue that more weight should be given to the larger N^{i) in a configuration since we are 
interested in the continuum physics and not in what happens in the stalk where N^i^i) is very 
small. We have tried various x^-fits with reasonable weights associated with the three- volumes 
N3{i). The kinetic term, which is the dominant term, is insensitive to any (reasonable) weight 
associated with N3{i). The potential term, which will be analyzed below, is more sensitive to 
the choice of the weight. However, the general power law dependence reported below is again 
unaffected by this choice. 
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values with the averages N^ii) measured directly. Similarly, we can use the mea- 
sured Mj's to reconstruct the second derivatives W'^N^^i)) and compare them to 



the form N-, 



-5/3/ 



coming from (19). 



The reconstruction of N-^^i) is illustrated in Fig. 4 for a variety of four- volumes 
iV4 and compared with the directly measured expectation values Nsi^i). It is 
seen that the reconstruction works very well and, most importantly, the coupling 
constant ki, which in this way is determined independently for each four- volume 
N4 really is independent of A^4 in the range of A^4's considered, as should be. 

We will now try to extract the potential [/"{N^li)) from the information con- 
tained in the matrix pp°*. The determination of (/"{N^li)) is not an easy task 
as can be understood from Fig. 5, which shows the measured coefficients Ui ex- 
tracted from the matrix pp°*, and which we consider somewhat remarkable. The 
interpolated curve makes an abrupt jump by two orders of magnitude going from 
the extended part of the universe (stretching over roughly 40 time steps) to the 
stalk. The occurrence of this jump is entirely dynamical, no distinction has ever 
been made by hand between stalk and bulk. 

There are at least two reasons for why it is difficult to determine the potential 
numerically. Firstly, the results are "contaminated" by the presence of the stalk. 
Since it is of cut-off size, its dynamics is dominated by fluctuations which likewise 
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Figure 5: Reconstructing the second derivative U"{N3{i)) from the coefficients 
Ui, for Ko = 2.2 and A = 0.6 and ivf '^^ = 160.000. 



are of cut-off size. They will take the form of short-time sub-dominant contri- 
butions in the correlator matrix C . Unfortunately, when we invert C to obtain 
the propagator P, the same excitations will correspond to the largest eigenvalues 
and give a very large contribution. Although the stalk contribution in the matrix 
C is located away from the bulk-diagonal, it can be seen from the appearance 
of the 1/A^^-term in eqs. (30) and (31) that after the projection orthogonal to 
the zero mode the contributions from the stalk will also affect the remainder of 
the geometry in the form of fluctuations around a small constant value. In de- 
riving Fig. 6 we have subtracted this constant value as best possible. However, 
the fluctuations of the stalk cannot be subtracted and only good statistics can 
eventually eliminate their effect on the behaviour of the extended part of the 
quantum universe. The second (and less serious) reason is that from a numerical 
point of view the potential term is always sub-dominant to the kinetic term for 
the individual spacetime histories in the path integral. For instance, consider the 
simple example of the harmonic oscillator. Its discretized action reads 



N 

Xj+i Xi \- 2 2 



s=j:At ( 



At 



(34) 



from which we deduce that the ratio between the kinetic and potential terms will 
be of order 1/At as At tends to zero. This reflects the well-known fact that the 
kinetic term will dominate and go to infinity in the limit as At -^ 0, with a typical 
path being nowhere differentiable. The same will be true when dealing with a 
more general action like (17) and its discretized version (19), where At scales like 
At ~ l/A'4 . Of course, a classical solution will behave differently: there the 
kinetic term will be comparable to the potential term. However, when extracting 
the potential term directly from the data, as we are doing, one is confronted with 
this issue. 

The range of the discrete three- volumes Ns{i) in the extended universe is from 
several thousand down to five, the kinematically allowed minimum. However, 
the behaviour for the very small values of N^li) near the edge of the extended 
universe is likely to be mixed in with discretization effects. In order to test 
whether one really has a N^ (i)-term in the action one should therefore only 
use values of iV3(i) somewhat larger than five. This has been done in Fig. 6, 
where we have converted the coefficients Ui from functions of the discrete time 
steps i into functions of the background spatial three- volume N^li) using the 
identification in (33) (the conversion factor can be read off the relevant curve in 
Fig. 4). It should be emphasized that Fig. 6 is based on data from the extended 
part of the spacetime only; the variation comes entirely from the central region 
between times -20 and 20 in Fig. 5, which explains why it has been numerically 
demanding to extract a good signal. The data presented in Fig. 6 were taken at a 
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Figure 6: The second derivative —U"{Ns) as measured for N^ 
Ko = 2.2 and A = 0.6. 



(4,1) 



160.000 and 



discrete volume A^4 = 160.000, and fit well the form N^ , corresponding to a 
potential ^2^3 • There is a very small residual constant term present in this fit, 
which presumably is due to the projection onto the space orthogonal to the zero 
mode, as already discussed earher. In view of the fact that its value is quite close 
to the noise level with our present statistics, we have simply chosen to ignore it 
in the remaining discussion. 

Apart from obtaining the correct power N^ for the potential for a given 
spacetime volume A^4, it is equally important that the coefficient in front of this 
term be independent of A'4. This seems to be the case as is shown in Fig. 7, 
where we have plotted the measured potentials in terms of reduced, dimensionless 
variables which make the comparison between measurements for different A^4's 
easier. - In summary, we conclude that the data allow us to reconstruct the 
action (19) with good precision. 

Let us emphasize a remarkable aspect of this result. Our starting point was 
the Regge action for CDT, as described in Sec. 2 above. However, the effective 
action we have generated dynamically by performing the nonperturbative sum 
over histories is only indirectly related to this "bare" action. Likewise, the cou- 
pling constant ki which appears in front of the effective action, and which we view 
as related to the gravitational coupling constant G by eq. (21) has no obvious 
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Figure 7: The dimensionless second derivative u = Nj U"{N^) plotted against 



N3/N, 



■3/4 



is the dimensionless spatial volume, for N, 



(4,1) 



40.000, 



i/~5/3, where u 

80.000 and 160.000, kq = 2.2 and A = 0.6. One expects a universal straight line 

near the origin (i.e. for large volumes) if the power law ^/(A^s) oc N^^^ is correct. 



direct relation to the "bare" coupling kq appearing in the Regge action (3) and 
in (7). Nevertheless the leading terms in the effective action for the scale factor 
are precisely the ones presented in (19). That a kinetic term with a second-order 
derivative appears as a leading term in an effective action is maybe less surprising, 
but it is remarkable and very encouraging for the entire CDT-quantization pro- 
gram that the kinetic term appears in precisely the correct combination with the 
factor A^3(z)^'^ needed to identify the leading terms with the corresponding terms 
in the Einstein-Hilbert action. In other words, only if these terms are present 
can we claim to have an effective field theory which has anything to do with the 
standard diffeomorphism-invariant gravitational theory in the continuum. This 
is neither automatic nor obvious, since our starting point involved both a dis- 
cretization and an explicit asymmetry between space and time, and since the 
nonperturbative interplay of the local geometric excitations we are summing over 
in the path integral is beyond our analytic control. Nevertheless, what we have 
found is that at least the leading terms in the effective action we have derived 
dynamically admit an interpretation as the standard Einstein term, thus passing 
a highly nontrivial consistency test. 
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Figure 8: Analyzing the quantum fluctuations of Fig. 1: diagonal entries F{t, t) 
of the universal scaling function F from (36), for iV] ' 
and 160.000. 
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5 Fluctuations around de Sitter space 

We have shown that the action (19) gives a very good description of the measured 
shape N^li) of the extended universe. Furthermore we have shown that by as- 
suming that the three- volume fluctuations around Nsii) are sufficiently small so 
that a quadratic approximation is valid, we can use the measured fluctuations to 
reconstruct the discretized version (19) of the minisuperspace action (17), where 
fci and /c2 are independent of the total four-volume A'^4 used in the simulations. 
This certainly provides strong evidence that both the minisuperspace descrip- 
tion of the dynamical behaviour of the (expectation value of the) three-volume, 
and the semiclassical quadratic truncation for the description of the quantum 
fluctuations in the three- volume are essentially correct. 

In the following we will test in more detail how well the actions (17) and (19) 
describe the data encoded in the covariance matrix C. The correlation function 
was defined in the previous section by 



C, 



nAi,i 



{6Nsit)6Ns{t')), SNsii) = N^it) - N^i{) 



(35) 



where we have included an additional subscript N4 to emphasize that A^4 is kept 
constant in a given simulation. The first observation extracted from the Monte 
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Carlo simulations is that under a change in the four- volume CN^(i,i') scales as'^ 

C^,{z, i') = N, f(i/nI'\ i'/nI") , (36) 

where F is a universal scaling function. This is illustrated by Fig. 8 for the 
rescaled version of the diagonal part Cj^^ (hi), corresponding precisely to the 
quantum fluctuations {{SN^i^i))'^)^/'^ of Fig. 1. While the height of the curve in 
Fig. 1 will grow as A^4 , the superimposed fluctuations will only grow as A^4 . 
We conclude that for fixed bare coupling constants the relative fluctuations will 
go to zero in the infinite-volume limit . 

From the way the factor y/N^ appears as an overall scale in eq. (20) it is 
clear that to the extent a quadratic expansion around the effective background 
geometry is valid one will have a scaling 

{SN3{t)6Ns{t')) = N!^'{Sn,{t,)6n,{t,)) = N,F{t„t,), (37) 

where tj = i/Nj . This implies that (36) provides additional evidence for the 
validity of the quadratic approximation and the fact that our choice of action 
(19), with ki independent of A^4 is indeed consistent. 

To demonstrate in detail that the full function F(t, t') and not only its diagonal 
part is described by the effective actions (17), (19), let us for convenience adopt a 
continuum language and compute its expected behaviour. Expanding (17) around 
the classical solution according to V^{t) = V^''{t)+x{t), the quadratic fluctuations 
are given by 

= M-\t,t'), (38) 

where T>x{s) is the normalized measure and the quadratic form M{t, t') is deter- 
mined by expanding the effective action S to second order in x{t), 

S{V,) = S{Vf) + -l^E.Jdtx{t)Hx{t). (39) 

In expression (39), H denotes the Hermitian operator 

^ _ d 1 d 4 

~ ~ dt cos3(t/5) dt ~ 52 cos5(t/5) ' ^ °^ 



^We stress again that the form (36) is only vahd in that part of the universe whose spatial 
extension is considerably larger than the minimal S'^ constructed from 5 tetrahedra. (The 
spatial volume of the stalk typically fluctuates between 5 and 15 tetrahedra.) 
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Figure 9: Comparing the two highest even eigenvector of the covariance matrix 
C(t, t') measured directly (gray curves) with the two lowest even eigenvectors of 
M~^(t,t'), calculated semiclassically (black curves). 

which must be diagonalized under the constraint that / dt^/g^ x{t) =0, since V4 
is kept constant. 

Let e*-"-'(t) be the eigenf unctions of the quadratic form given by (39) with the 
volume constraint enforced^, ordered according to increasing eigenvalues A„. As 
we will discuss shortly, the lowest eigenvalue is Ai = 0, associated with transla- 
tional invariance in time direction, and should be left out when we invert M{t, t'), 
because we precisely fix the centre of volume when making our measurements. 
Its dynamics is therefore not accounted for in the correlator C{t,t'). 

If this cosmological continuum model were to give the correct description of 



*^One simple way to find the eigenvalues and eigenfunctions approximately, including the con- 
straint, is to discrctize the differential operator, imposing that the (discretized) eigenfunctions 

vanish at the boundaries t = ±_B7r/2 and finally adding the constraint as a term H J dtx{t) J 
to the action, where the coefficient ^ is taken large. The differential operator then becomes 
an ordinary matrix and eigenvalues and eigenvectors can be found numerically. Stability with 
respect to subdivision and choice of ^ is easily checked. 
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the computer-generated universe, the matrix 



M-\t,t') = J2 , ^ • (41) 



n=2 



should be proportional to the measured correlator C{t,t'). Fig. 9 shows the 
eigenfunctions e*-^-'(t) and e^^\t) (with two and four zeros respectively), calculated 
from H with the constraint J dt^/g^ x{t) = imposed. Simultaneously we show 
the corresponding eigenfunctions calculated from the data, i.e. from the matrix 
C{t,t'), which correspond to the (normalizable) eigenfunctions with the highest 
and third-highest eigenvalues. The agreement is very good, in particular when 
taking into consideration that no parameter has been adjusted in the action 
(we simply take B = SqNJ At in (14) and (39), which gives B = 14.47aj for 
Ni = 362.000). 

The reader may wonder why the first eigenfunction exhibited has two zeros. 
As one would expect, the ground state eigenfunction e*^°-*(t) of the Hamiltonian 
(40), corresponding to the lowest eigenvalue, has no zeros, but it does not satisfy 
the volume constraint J dty/g^ x{t) = 0. The eigenfunction e^^\t) of H with 
next-lowest eigenvalue has one zero and is given by the simple analytic function 

where c is a constant. One realizes immediately that e*^^-* is the translational zero 
mode of the classical solution V^\t) (oc cos^ t/B). Since the action is invariant 
under time translations we have 

S{Vf{t + At)) = S{Vf{t)), (43) 

and since V^\t) is a solution to the classical equations of motion we find to second 
order (using the definition (42)) 

S{Vf{t + At)) = S{Vf{t)) + 4^^ / dt e^'\t)He^'\t), (44) 

consistent with e^^\t) having eigenvalue zero. 

It is clear from Fig. 9 that some of the eigenfunctions of H (with the volume 
constraint imposed) agree very well with the measured eigenfunctions. All even 
eigenfunctions (those symmetric with respect to reflection about the symmetry 
axis located at the centre of volume) turn out to agree very well. The odd 
eigenfunctions of H agree less well with the eigenfunctions calculated from the 
measured C{t, t'). The reason seems to be that we have not managed to eliminate 
the motion of the centre of volume completely from our measurements. As already 
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mentioned above, there is an inherent ambiguity in fixing the centre of volume, 
which turns out to be sufficient to reintroduce the zero mode in the data. Suppose 
we had by mistake misplaced the centre of volume by a small distance At. This 
would introduce a modification 

A,/3 = ^ At (45) 

dt ^ ^ 

proportional to the zero mode of the potential V^''{t). It follows that the zero 
mode can re-enter whenever we have an ambiguity in the position of the centre 
of volume. In fact, we have found that the first odd eigenfunction extracted 
from the data can be perfectly described by a linear combination of e'^^\t) and 
e'^^\t). It may be surprising at first that an ambiguity of one lattice spacing 
can introduce a significant mixing. However, if we translate AV3 from eq. (45) 
to "discretized" dimensionless units using V3(i) ~ N^' cos{i/Nj ), we find that 
AV3 ~ \/iV4, which because of {{SN^^i))^) ~ N4 is of the same order of magnitude 
as the fluctuations themselves. In our case, this apparently does affect the odd 
eigenfunctions. 

One can also compare the data and the matrix M~^{t, t') calculated from (41) 
directly. This is illustrated in Fig. 10, where we have restricted ourselves to data 
from inside the extended part of the universe. We imitate the construction (41) 
for M~^, using the data to calculate the eigenfunctions, rather than H. One 
could also have used C(t, t') directly, but the use of the eigenfunctions makes it 
somewhat easier to perform the restriction to the bulk. The agreement is again 
good (better than 15% at any point on the plot), although less spectacular than 
in Fig. 9 because of the contribution of the odd eigenfunctions to the data. 

6 The geometry of spatial three-spheres 

We have shown above that our data for the spatial three-volumes have a natural 
interpretation as coming from the slicing of a four-sphere with standard geometry 
(the "round" four-sphere), with relatively small quantum ffuctuations superim- 
posed. It is natural to ask to what extent the spatial three-spheres themselves 
can be assigned the standard geometry of a "round" three-sphere, again with 
relatively small quantum ffuctuations superimposed. We have already provided 
evidence that the Hausdorff dimension of the spatial slices is three [6, 4]. How- 
ever, the Hausdorff dimension is a very coarse measure of geometry, and even 
very fractal structures can have Hausdorff dimension three ^. 



^To illustrate the point, the Hausdorff dimension of the complex plane (with standard geom- 
etry) is of course equal to two, but the same is true for the highly fractal structure of so-called 
branched polymers or planar trees embedded in the plane. 
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Figure 10: Comparing data for the extended part of the universe: measured 
C{t,t') (above) versus M~^{t,t') obtained from analytical calculation (below). 
The agreement is good, and would have been even better had we included only 
the even modes. 

We have analyzed the geometry of the spatial three-spheres as follows. Each 
spatial slice at integer proper time i is a triangulation, consisting of a certain 
number N3 of tetrahedra, glued together pairwise such that the resulting topology 
is that of a three-sphere. We now choose an arbitrary tetrahedron as the origin of 
measurements and subsequently decompose the S^ into (thick) shells of tetrahedra 
characterized by their distance r from this origin, where the distance r is defined 
as the minimal number of tetrahedra one has to cross when moving from the 
shell to the origin via neighbouring tetrahedra. We call the number of tetrahedra 
in the shell at distance r the area A{r, N3) of the shell. In order to compute 
the expectation value of this quantity, we have to repeat the measurements in a 
way that averages over different triangulations of S^, over different spatial slices 
within the S'^'s and over different locations of the point of origin within those 
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X = rlNl'^ii) 

Figure 11: Testing relation (46) for the bare coupling constants kq = 2.2 and 
A = 0.6, at four-volume AQ = 160.000. Data have been collected for spatial 
slices at various distances close to the centre of volume. 

slices. In this manner we can test whether {A{r,N^)) behaves like a regular 
three-sphere (with only small fluctuations superimposed), with r viewed as the 
geodesic distance. If this was the case, one would expect a functional dependence 
of the form 

(46) 



{A{r,Ns))^N^^' 



sm 



cN. 



1/3 



with c a constant. 

Fig. 11 summarizes the results of our measurements. Since we are not in- 
terested in very small iVa's where no continuum scaling is expected, we have re- 
stricted ourselves to spatial shces close to the centre of volume as defined above, 
where N^ is largest. The first thing to note about Fig. 11 is that the data from 
spatial slices at different distances from the centre of volume fall to good accuracy 
on a common, universal curve. Next, we observe that relation (46) is reasonably 
well satisfied, except for the measurements at large radii r, which exhibit a tail 
not described by formula (46). This signals the presence of large fluctuations in 
the geometry (the shape) of the spatial slices to the effect that we cannot simply 
view them - in the sense of expectation values - as classical spheres of constant 
positive curvature with fixed radius proportional to N'^^^, superimposed by small 
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quantum fluctuations. In fact there is already evidence that the geometry, when 
defined with respect to the geodesic distance r, has certain fractal properties [4]. 
This can be substantiated and quantified by measuring the topology of a typical 
spherical shell at a distance r from a chosen origin in more detail. At sufficiently 
large radius r one finds that the topology is no longer that of a single two-sphere, 
but branches out into a number of disconnected pieces, most likely by effectively 
creating a number of spatial "baby universes" . It is well known how to study the 
distribution of such baby universes [25, 26], and we believe that these methods 
will yield a quantitative description of the observed slower fall-off for large r. De- 
tails of this picture, including a study of the temporal dynamics of such spatial 
baby universes, will be published elsewhere. 

7 The size of the universe and the flow of G 

Let us now return to equation (21), 

G = ^^^, (47) 

which relates the parameter ki extracted from the Monte Carlo simulations to 
Newton's constant in units of the cut-off a, G/a^. For the bare coupling con- 
stants (kq. A) = (2.2, 0.6) we have high-statistics measurements for N^^ ranging 
from 45.500 to 362.000 four-simplices (equivalently, AQ ' ranging from 20.000 to 
160.000 four-simplices). The choice of A determines the asymmetry parameter a, 
and the choice of (kq, A) determines the ratio ^ between A^] ' and AQ . This in 
turn determines the "effective" four-volume C4 of an average four-simplex, which 
also appears in (47). The number sq in (47) is determined directly from the time 
extension Tuniv of the extended universe according to 

Tuniv = 7r5o(A^i'''^)'^'. (48) 

Finally, from our measurements we have determined ki = 0.038. Taking every- 
thing together according to (47), we obtain G ~ 0.23a^, or ipi ^ 0.48a, where 
ipi = vG is the Planck length. 

From the identification of the volume of the four-sphere, V4 = Sn'^R^/S = 
C^N^ ' a^, we obtain that R = 3.1a. In other words, the linear size irR of the 
quantum de Sitter universes studied here lies in the range of 12-21 Planck lengths 
for A4 in the range mentioned above and for the hare coupling constants chosen 
as (fi;o,A) = (2.2, 0.6). 1° 



^° Small deviations from the corresponding numbers quoted in [8] have their origin in the 
more careful (and correct) treatment of the various four-volumes N^, AQ ' and TV^ ' in the 
present work. 
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Our dynamically generated universes are therefore not very big, and the quan- 
tum fluctuations around their average shape are large as is apparent from Fig. 
1. It is rather surprising that the semiclassical minisuperspace formulation is 
applicable for universes of such a small size, a fact that should be welcome news 
to anyone performing semiclassical calculations to describe the behaviour of the 
early universe. However, in a certain sense our lattices are still coarse compared 
to the Planck scale ipi because the Planck length is roughly half a lattice spacing. 
If we are after a theory of quantum gravity valid on all scales, we are in particu- 
lar interested in uncovering phenomena associated with Planck-scale physics. In 
order to collect data free from unphysical short-distance lattice artifacts at this 
scale, we would ideally like to work with a lattice spacing much smaller than the 
Planck length, while still being able to set by hand the physical volume of the 
universe studied on the computer. 

The way to achieve this, under the assumption that the coupling constant G 
of formula (47) is indeed a true measure of the gravitational coupling constant, 
is as follows. We are free to vary the discrete four-volume A^4 and the bare 
couphng constants (kq, A) of the Regge action (see [4] for further details on the 
latter). Assuming for the moment that the semiclassical minisuperspace action 
is valid, the effective coupling constant ki in front of it will be a function of the 
bare coupling constants (kq, A), and can in principle be determined as described 
above for the case (kq, A) = (2.2, 0.6). If we adjusted the bare couphng constants 
such that in the limit as A^4 -^ oo both 

Vi^N^a^ and G ^a^/h{Ko,A) (49) 

remained constant (i.e. /ci(fi;o. A) ~ I/a/A^), we would eventually reach a region 
where the Planck length was significantly smaller than the lattice spacing a, in 
which event the lattice could be used to approximate spacetime structures of 
Planckian size and we could initiate a genuine study of the sub-Planckian regime. 
Since we have no control over the effective coupling constant ki, the first obvious 
question which arises is whether we can at all adjust the bare coupling constants 
in such a way that at large scales we still see a four- dimensional universe, with 
ki going to zero at the same time. The answer seems to be in the affirmative, as 
we will go on to explain. 

Fig. 12 shows the results of extracting ki for a range of bare coupling constants 
for which we still observe an extended universe. In the top figure A = 0.6 is kept 
constant while kq is varied. For kq sufficiently large we eventually reach a point 
where a phase transition takes place (the point in the square in the bottom right- 
hand corner is the measurement closest to the transition we have looked at). For 
even larger values of kq, beyond this transition, the universe disintegrates into a 
number of small universes, in a CDT-analogue of the branched-polymer phase of 
Euclidean quantum gravity. The plot shows that the effective coupling constant 
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Figure 12: The measured effective coupling constant ki as function of tlie bare kq 
(top, A = 0.6 fixed) and the asymmetry A (bottom, kq = 2.2 fixed). The marked 
point near the middle of the data points sampled is the point (kq. A) = (2.2, 0.6) 
where most measurements in the remainder of the paper were taken. The other 
marked points are those closest to the two phase transitions, to the "branched- 
polymer phase" (top), and the "crumpled phase" (bottom). 

ki becomes smaller and possibly goes to zero as the phase transition point is 
approached, although our current data do not yet allow us to conclude that ki 
does indeed vanish at the transition point. 

Conversely, the bottom figure of Fig. 12 shows the effect of varying A, while 
keeping kq = 2.2 fixed. As A is decreased towards 0, we eventually hit another 
phase transition, separating the physical phase of extended universes from the 
CDT-equivalent of the crumpled phase of Euclidean quantum gravity, where the 
entire universe will be concentrated within a few time steps, as already mentioned 
in Sec. 3 above. (The point closest to the transition where we have taken mea- 
surements is the one in the bottom left-hand corner.) Also when approaching 
this phase transition the effective coupling constant ki goes to 0, leading to the 
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tentative conclusion that fci — > along the entire phase boundary. 

However, to extract the coupling constant G from (47) we not only have 
to take into account the change in /ci, but also that in sq (the width of the 
distribution N^{i)) and in the effective four- volume C4 as a function of the bare 
coupling constants. Combining these changes, we arrive at a slightly different 
picture. Approaching the boundary where spacetime collapses in time direction 
(by lowering A), the gravitational coupling constant G decreases, despite the fact 
that 1/ki increases. This is a consequence of §0 decreasing considerably, as can 
be seen from Fig. 3. On the other hand, when (by increasing kq) we approach 
the region where the universe breaks up into several independent components, 
the effective gravitational coupling constant G increases, more or less like l/Zci, 
where the behaviour of ki is shown in Fig. 12 (top). This implies that the Planck 
length ^pi = vG increases from approximately 0.48a to 0.83a when kq changes 
from 2.2 to 3.6. Most likely we can make it even bigger in terms of Planck units 
by moving closer to the phase boundary. 

On the basis of these arguments, it seems likely that the nonperturbative 
CDT-formulation of quantum gravity does allow us to penetrate into the sub- 
Planckian regime and probe the physics there explicitly. Work in this direction 
is currently ongoing. One interesting issue under investigation is whether and 
to what extent the simple minisuperspace description remains valid as we go to 
shorter scales. We have already seen deviations from classicality at short scales 
when measuring the spectral dimension [20, 4], and one would expect them to 
be related to additional terms in the effective action (17) and/or a nontrivial 
scaling behaviour of ki. This raises the interesting possibility of being able to 
test explicitly the scaling violations of G predicted by renormalization group 
methods in the context of asymptotic safety [2]. 

8 Discussion 

The CDT model of quantum gravity is extremely simple. It is the path integral 
over the class of causal geometries with a global time foliation. In order to perform 
the summation explicitly, we introduce a grid of piecewise linear geometries, much 
in the same way as when defining the path integral in quantum mechanics. Next, 
we rotate each of these geometries to Euclidean signature and use as bare action 
the Einstein-Hilbert action^^ in Regge form. That is all. 

The resulting superposition exhibits a nontrivial scaling behaviour as func- 
tion of the four-volume, and we observe the appearance of a well-defined average 
geometry, that of de Sitter space, the maximally symmetric solution to the clas- 



^^Of course, the full, effective action, including measure contributions, will contain all higher- 
derivative terms. 
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sical Einstein equations in the presence of a positive cosmo logical constant. We 
are definitely in a quantum regime, since the fluctuations of the three-volume 
around de Sitter space are sizeable, as can be seen in Fig. 1. Both the average 
geometry and the quantum fluctuations are well described in terms of the mini- 
superspace action (17). A key feature to appreciate is that, unlike in standard 
(quantum-) cosmological treatments, this description is the outcome of a nonper- 
turbative evaluation of the full path integral, with everything but the scale factor 
(equivalently, V3(t)) summed over. Measuring the correlations of the quantum 
fluctuations in the computer simulations for a particular choice of bare coupling 
constants enabled us to determine the continuum gravitational coupling constant 
G as G ~ 0.42a^, thereby introducing an absolute physical length scale into the 
dimensionless lattice setting. Within measuring accuracy, our de Sitter universes 
(with volumes lying in the range of 6.000-47.000 ipi) are seen to behave perfectly 
semiclassically with regard to their large-scale properties. 

We have also indicated how we may be able to penetrate into the sub-Planckian 
regime by suitably changing the bare coupling constants. By "sub-Planckian 
regime" we mean that the lattice spacing a is (much) smaller than the Planck 
length. While we have not yet analyzed this region in detail, we expect to eventu- 
ally observe a breakdown of the semiclassical approximation. This will hopefully 
allow us to make contact with attempts to use renormalization group techniques 
in the continuum and the concept of asymptotic safety to study scaling violations 
in quantum gravity [2]. 

On the basis of the results presented here, two major issues suggest themselves 
for further research. First, we need to establish the relation of our effective 
gravitational coupling constant G with a more conventional gravitational coupling 
constant, defined directly in terms of coupling matter to gravity. In the present 
work, we have defined G as the coupling constant in front of the effective action, 
but it would be desirable to verify directly that a gravitational coupling defined 
via the coupling to matter agrees with our G. In principle it is easy to couple 
matter to our model, but it is less straightforward to define in a simple way a set- 
up for extracting the semiclassical effect of gravity on the matter sector. Attempts 
in this direction were already undertaken in the "old" Euclidean approach [22, 23], 
and it is possible that similar ideas can be used in CDT quantum gravity. Work 
on this is in progress. 

The second issue concerns the precise nature of the "continuum limit" . Recall 
our discussion in the Introduction about this in a conventional lattice-theoretic 
setting. The continuum limit is usually linked to a divergent correlation length 
at a critical point. It is unclear whether such a scenario is realized in our case. In 
general, it is rather unclear how one could define at all the concept of a divergent 
length related to correlators in quantum gravity, since one is integrating over all 
geometries, and it is the geometries which dynamically give rise to the notion of 
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"length". 

This has been studied in detail in two-dimensional (Euclidean) quantum grav- 
ity coupled to matter with central charge c < 1 [24] . It led to the conclusion that 
one could associate the critical behaviour of the matter fields (i.e. approaching 
the critical point of the Ising model) with a divergent correlation length, al- 
though the matter correlators themselves had to be defined as non-local objects 
due to the requirement of diffeomorphism invariance. On the other hand, the 
two-dimensional studies do not give us a clue of how to treat the gravitational 
sector itself, since they do not possess gravitational field-theoretic degrees of free- 
dom. What happens in the two-dimensional lattice models which can be solved 
analytically is that the only fine-tuning needed to approach the continuum limit 
is an additive renormalization of the cosmological constant (for fixed matter cou- 
plings). Thus, fixing the two-dimensional spacetime volume A''2 (the number of 
triangles), such that the cosmological constant plays no role, there are no further 
coupling constants to adjust and the continuum limit is automatically obtained 
by the assignment V2 = N2a'^ and taking A^2 ^ 00. This situation can also occur 
in special circumstances in ordinary lattice field theory. A term like 

^Ci(0i+i - 00^ + C2(0i+1 + 0i-l - 20i)^ (50) 

i 

(or a higher- dimensional generalization) will also go to the continuum free field 
theory simply by increasing the lattice size and using the identification Vd = L'^a'^ 
{L denoting the linear size of the lattice in lattice units), the higher-derivative 
term being sub-dominant in the limit. It is not obvious that in quantum grav- 
ity one can obtain a continuum quantum field theory without fine-tuning in a 
similar way, because the action in this case is multiplied by a dimensionful cou- 
pling constant. Nevertheless, it is certainly remarkable that the infrared limit 
of our effective action apparently reproduces - within the cosmological setting - 
the Einstein-Hilbert action, which is the unique diffeomorphism-invariant gener- 
alization of the ordinary kinetic term, containing at most second derivatives of 
the metric. A major question is whether and how far our theory can be pushed 
towards an ultraviolet limit. We have indicated how to obtain such a limit by 
varying the bare coupling constants of the theory, but the investigation of the 
limit a ^ with fixed G has only just begun. 
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